function removeD = removeDigonal(M)
M(logical(eye(size(M)))) = 0;
b=triu(M);
b=b(:,2:end);
c=tril(M);
c=c(:,1:end-1);
removeD = b+c;
